Superconducting proximity effect in graphene under inhomogeneous strain 
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The interplay between quantum Hall states and Cooper pairs is usually hindered by the suppres- 
sion of the superconducting state due to the strong magnetic fields needed to observe the quantum 
Hall effect. From this point of view graphene is special since it allows the creation of strong pseudo- 
magnetic fields due to strain. We show that in a Josephson junction made of strained graphene, 
Cooper pairs will diffuse into the strained region. The pair correlation function will be sub-lattice 
polarized due to the polarization of the local density of states in the zero pseudo-Landau level. We 
uncover two regimes; 1) one in which the cyclotron radius is larger than the junction length in which 
case the supercurrent will be enhanced, and 2) the long junction regime where the supercurrent is 
strongly suppressed because the junction becomes an insulator. In the latter case quantized Hall 
states form and Andreev scattering at the normal/superconducting interface will induce edge states. 
Our numerical calculation has become possible due to an extension of the Chebyshev Bogoliubov-de 
Gennes method to computations on video cards (CPUs). 

PACS numbers: 73.22. Pr, 74.45. -|-c,73.43.-f 



Graphene has eome in reeent years to the forefront 
of condensed matter research not only due to its pecu- 
liar electronic properties but also due to its technological 
potential^ The special arrangement of carbon atoms 
in a honeycomb lattice has major consequences for its 
electronic properties. The particular way of coupling be- 
tween lattice deformations and the electronic states re- 
sults in remarkable properties. Any strain (either intrin- 
sically due to phonons'^ or ripples^ or extrinsically due to 
applied stress^"^° couples to the electronic degrees of free- 
dom as a gauge field. The induced field does not break 
time reversal symmetry, because it has opposite sign in 
the K and K' Dirac cones. For this reason the field is 
referred to as a pseudo-magnetic field. Electrons with 
momentum pertaining to different valleys will feel the 
effect of a reversed magnetic field by being deflected in 
opposite directions'^. More interestingly, it was recently 
shown^'^ that when the pseudo-magnetic field is slowly 
varying, the electronic spectrum shows the appearance 
of pseudo-Landau levels. Another special ingredient in 
our study is the existence of superconducting correla- 
tions in the graphene layer. While intrinsic supercon- 
ductivity was not observed experimentally and only pre- 
dicted theoretically in very specific conditions, e.g. the 
presence of high doping by adatoms'^, superconducting 
correlations can be induced by proximity to a supercon- 
ducting contact"'^'^"'^. Putting together these two effects 
one can achieve the coexistence of pseudo-quantum Hall 
states and superconducting correlations. Of special inter- 
est are the edge states formed at the interface between 
normal and superconducting regions since for energies be- 
low the superconducting gap quasiparticles will undergo 
Andreev scattering^''. The coexistence between super- 
conductivity and quantum Hall states was observed in 
superconductor/2DEG structures in magnetic fields. In 
order to observe quantized conductances the magnetic 
length has to be small when compared to the scattering 
length but still large when compared to the supercon- 
ducting coherence length'^"^°. The coexistence between 
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FIG. 1. (Color online) Strained graphene Josephson junc- 
tion obtained from a rectangle by applying linearly varying 
forces along the edges. The grayscale (color) code of the inter- 
atomic links is given by the ratio tij/7o and shows that the 
links along x are stretched while the ones along y are com- 
pressed. Depending on the strength of the pseudo-magnetic 
field the cyclotron radius is such that it can take the quasi- 
particle across the junction or generate an edge state. 



antagonistic states can be achieved in graphene because 
the pseudo-magnetic field couples to the electronic or- 
bital degrees of freedom in the graphene layer while leav- 
ing the superconducting contacts unaffected. One can 
envision the generation of edge states generated by large 
pseudo-magnetic fields (even on the order of hundreds 
of Tesla) but which undergo Andreev reflections at the 
superconducting/graphene interface. From this point of 
view, graphene is unique. In this paper we describe a 
strained graphene Josephson junction and show how the 
proximity effect can be tuned by inhomogeneously strain- 
ing the junction and how the Josephson current could be 
carried either across the junction or by edge states. 



We consider a Josephson junction made of strained 
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graphene. Previous theoretical descriptions of graphene 
Joshepson junctions used both the low energy Dirac 
formulation^^ and tight- binding Hamiltonians^^ . Our 
model follows closely Ref. (22) which consider that 
graphene under the contacts has an intrinsic pairing po- 
tential, then the Bogoliubov-de Gennes equations are 
solved self-consistently. To solve the problem we take 
instead another route by considering three dimensional 
superconducting contacts. Since the inverse proximity 
effect of the graphene layer on the contacts will be neg- 
ligible and since there are no intrinsic superconducting 
correlations in graphene, a non self-consistent approach is 
reasonable. Our approach is also based on tight-binding 
Hamiltonians, therefore it will give similar results for the 
unstrained junction and not shine light in the conflict- 
ing results of Refs. (21) and (22). The system shown in 
Fig. (1) has contacts made of multilayer (10 layers) AA 
stacked graphite which is strongly coupled to the two- 
dimensional graphene sheet. Also from an experimental 
point of view, in order to have a sizable proximity ef- 
fect the contact resistance between the contacts and the 
graphene layer has to be very low. For the purpose of our 
calculation, e.g. the superconducting proximity effect, 
the specifics of the contact band structure and coupling 
to the graphene sheet are not important. A detailed cal- 
culation which takes into account the coupling between a 
metallic gate and graphene could in principle be done but 
the quantitative description of the contacts will not add 
much in terms of qualitative description of the system. 
Such strained junctions could be experimentally achieved 
by rotating the superconducting contacts with respect to 
the symmetry axis of the junction. 

The electronic properties are described by a tight- 
binding Hamiltonian for the tt carbon orbitals. The 
minimal Hamiltonian needed to describe superconduct- 
ing correlations has the following Nambu spinor form: 




where T-Lij is a 2 x 2 matrix: 

H.^(^5^^')^. + (;^-|)(l-%). (2) 

where the sum is over nearest and next nearest 

neighbors. The superconducting order parameter, A^, is 
of s-wave spin-singlet type and is non-zero only in the 
contact region above the graphene sheet. The strain 
information is included in the modified hopping am- 
plitudes tij according to the empirical relation tij = 

7o exp^''^''^~ ^\ where 70 = 3eV and oq is the equi- 
librium inter-carbon distance^'^. This also gives a good 
approximation for the next-nearest neighbor hopping am- 
plitude. In order to describe the inhomogeneous strain 
we consider the deformed rectangular geometry shown in 
Fig. (1) which can be obtained from a rectangle by ap- 
plying a linearly varying force = {—const * y, 0) on 



the left side and Fr = {const * y, 0) on the right side^. 
We obtain the strain tensor by numerically solving the 
elasticity problem for an isotropic elastic sheet with the 
elastic properties of an ideal flat graphene sheet. 

Since we consider three dimensional contacts and the 
strain is inhomogeneous it is not possible to use any sym- 
metries in the calculations of the electronic properties. 
Therefore the system size needed to describe this prob- 
lem becomes prohibitively large for an exact diagonal- 
ization of the Hamiltonian. Instead we will numerically 
obtain an approximation of the Gorkov Green's function 
by using the Chebyshev Bogoliubov de Gennes (CBdG) 
method^-*. The 2x2 Green's function is defined as: 




where G{uj + irj) ~ [Ld+i7]—7{\~^, \vac) is the vacuum and 
the diagonal(normal) and off-diagonal(anomalous) com- 
ponents can be expressed as: 

Gl}{u;) = (c.t|G(..)|4^) (4) 

Glf{u:)^{cl\G{u:)\c]^r. (5) 

First a scaling of the excitation energies is performed, 
e.g. H = {H — lb) /a, uj = {oj — b)/a where a = {E^ax — 
E„itn)/ (2 - 77) and b = {Ej^ax + Ermn)/2, where 77 > is 
a small number. Following Refs. [24 and 25], the Green's 
function's components can be expressed as an expansion 
written in terms of Chebyshev polynomials: 

Gl^^'^iu) = ^!^^^«ii(i2)(,,,)e— cosM (6) 

^ n—0 

where the coefficients al^'"^^\i,j) can be obtained by 
an iterative procedure involving repeated applications 
of the Hamiltonian on iterative vectors. The physi- 
cal properties that can be straightforwardly extracted 
from the Green's functions are the local density of states 
(LDOS) Ni{uj) = -^Im[G}l{uj)] and the spin-singlet 
superconducting pair correlation function < Ci-^Cii >= 

^^XG\f{E){l-2f{E))dE. 

Even though it is feasible to use the CBdG method 
on parallel computer clusters even for a large number of 
atoms^"* (even up to one million) we have instead devel- 
oped codes which take advantage of the parallel nature 
of Graphics Processing Units(GPU). Significant speedups 
on the order of 1000 can be achieved for calculating elec- 
tronic properties of systems described by tight-binding 
Hamiltonians; the system size is on the order of hun- 
dreds of thousands of atoms. The calculations presented 
here have been performed on a system containing three 
Nvidia Geforce GTX580 which ran in parallel. 



In Fig. (2a) we show the pair correlation function as a 
function of position for various strain configurations for 
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FIG. 2. (Color online) a) Absolute value of the pair corre- 
lation function along a line, y — Ly/2, as a function of po- 
sition for various strain configurations for /i — O.2470 (the 
zero pseudo-Landau level is at the Fermi level). The inset 
shows the pair correlation function for fi = 0. b) Local den- 
sity of states along the same line as a function of position for 
/i = 0.2470. The inset shows the LDOS for fj, = 0. 




FIG. 3. (Color online) Pair correlation at a; = Lx/2, y = Ly/2 
as a function of for different strain configurations. For the 
A sub-lattice two regimes are found depending on the ratio 
Tc/Lx. Inset: the same but now for the B sub-lattice where 
the pair correlation decays exponentially with Lx- 

in strained Josephson junctions are expected; flow will 
only occur in one sub-lattice (A). The inset of Fig. (2b) 
shows the LDOS for /i = and since the Fermi level is 
away from the zero pseudo-Landau level there is no sym- 
metry breaking and the proximity effect is a conventional 
one. 



jjL = 0.2470. The plot is taken along a line y = Ly/2. 
The maximum strain, max{uxx) = £maxi gives the cur- 
vature of the deformed junction and is a direct measure 
of the pseudo-magnetic field ""^^ . As expected, under the 
superconducting contact the pair correlation is constant 
while in the junction it decreases exponentially. No- 
tice that as the strain increases the results for the two 
sub-lattices start to differ. While for one sub-lattice the 
pair correlation is enhanced when compared to the un- 
strained junction, in the other sub-lattice it is strongly 
suppressed. The chemical potential is chosen such that 
the zero pseudo-Landau level is pinned at the Fermi level 
(since we also consider next nearest neighbors it is shifted 
away from ji = 0). The inset shows the same quantity 
but for ^ = and shows no significant difference on the 
two sub-lattices. In order to understand this peculiar be- 
havior we also plot in Fig. (2b) the LDOS for the same 
configurations. We see a breaking of the sub-lattice sym- 
metry which shows that the LDOS is large in one sub- 
lattice and suppressed in the other. The symmetry break- 
ing of the zcro-th pseudo-Landau level can be observed 
even in a strained graphene sheet with no superconduct- 
ing contacts and is a direct consequence of the opening 
of a gap. It can be induced by breaking either the time 
reversal symmetry or the sub-lattice symmetry and since 
the pseudo-magnetic fields have no net flux (they have al- 
ternating sign in the two Dirac cones) the opening of the 
gap can be achieved by breaking the symmetry of the zero 
pseudo-Landau level. Changes in the super-current flow 



The dependence of the pair correlation at the mid 
point {Lx/2, Ly/2) for fi = O.2470 as a function of the 
length of the junction is shown in Fig. (3) where two 
regimes are found. First, when the junction is shorter 
than twice the effective cyclotron radius, Tc oc ^y/ tmax-, 
quasi-particlcs with energy lower than the superconduct- 
ing gap are bound and after undergoing Andreev scat- 
tering at the normal/superconducting interface they will 
scatter between the contacts. In this regime the pair 
correlation in the A sub-lattice is enhanced, because the 
LDOS is itself enhanced by the pseudo-quantum Hall ef- 
fect. The diffusion length of Cooper pairs is larger there- 
fore enhancing the pair correlations at Lx/2. Second, 
when the junction is longer than twice the effective cy- 
clotron radius quantized states will form in the junction 
which becomes insulating. Andreev scattering will only 
affect edge states near the normal/superconducting in- 
terface. In this regime the diffusion length is strongly 
suppressed and the pair correlation at Lx/2 shows an 
exponential decay. The inset of Fig. (3) shows that in 
the B sub-lattice the diffusion length of the Cooper pairs 
decreases monotonically with the strain and the pair cor- 
relation at Lx/2 decays exponentially with L^ faster than 
the pair correlation in the A sub-lattice due to the ab- 
sence of states localized in the B sub-lattice. 

These two regimes are also present when considering 
the formation of quantum Hall states in the presence of 
pseudo-magnetic fields. This is shown in Fig. (4) where 
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FIG. 4. (Color online) LDOS at the Fermi level for the A 
sub-lattice across the junction aX y — Ly/2 for different Lx- 

we plot the LDOS at the Fermi level for the A sub-lattice 
for the same strain configuration, emax = 0.2, but for 
different lengths of the junction. For short junctions the 
quasi-particles scatter from both contacts and the forma- 
tion of quantized orbits is hindered. For long junctions 
we observe a plateau in the LDOS signaling the appear- 
ance of quantized states in the center of the junctions. 
In this case the Andreev scattering affects only a narrow 
region near the normal/superconducting interface. 



To further see how the pseudo-Landau levels disperse 
near the interface we plot in Fig. (5) the LDOS across 
the junction for fixed y = Ly/2 and atoms in the A 
sub-lattice for three strain configurations. When the 
pseudo-magnetic field is strong {emax — 0.3) we ob- 
serve the formation of pseudo-Landau levels (LL) (up to 
ri = 5) which are localized in the junction. They disperse 
only very close to the normal/superconducting junctions 
{Lx/ao = 75 and Lx/ao = 175) and even show quan- 
tized edge states due to Andreev reflections, i.e. these 
are Andreev Bound States (ABS). In the B sub- lattice 
the LDOS is similar except at the Fermi level where the 
zero pseudo-Landau level is missing. As the strain is re- 
duced to emax = 0.2 and emax ~ 0.15 the pseudo-Landau 
levels become dispersive away from the interfaces and in 
the latter case quantized states due to Andreev reflec- 
tions can be observed throughout the whole junction. 



In conclusion we have shown that by using an efficient 
numerical method (CBdG) which is extended to run on 
CPUs with the benefit of large speedups, simulations of 
strained grap Irene junctions with a large number or atoms 
are possible. We showed for the first time that strained 
graphene proves to be a playground for testing the ex- 
otic coexistence between strong pseudo-magnetic fields 



and superconducting correlations. This is due to the fact 
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FIG. 5. (Color online) LDOS across the junction at y = Ly/2 
for a) emax ~ 0.3 and b) emax ~ 0.1. The value of the order 
parameter in the contacts is A = O.270 and the chemical 
potential is /i = O.2470. As a guide to the eye we show the 
position of the LL(ABS) by horyzontal black(red) lines and 
the superconducting gap (A). 



that the pseudo-magnetic fields only couple to the or- 
bital motion of electrons in the graphene sheet and not 
to their spin. Any superconducting contact will be un- 
affected by these strong pseudo-magnetic fields. By con- 
sidering a Josephson junction made of strained graphene 
we showed that due to the proximity effect Cooper pairs 
can leak into the strained region. In addition, due to 
the peculiar nature of the pseudo-Landau levels the pair 
correlation function has a broken sub-lattice symmetry, 
being enhanced in one sub-lattice and suppressed in the 
other. We have also uncovered two regimes depending on 
the ratio between the cyclotron radius and the junction 
length; in one regime quasi-particles will bounce between 
the superconducting contacts in which case the Joseph- 
son current would be enhanced while in the other regime 
quasi-particles scatter only from one contact thus cre- 
ating an Andreev bound edge state in which case the 
Josephson current would be strongly suppressed across 
the junction. In this case the supercurrent can only 
propagate through edge states. From an experimental 
point of view, the particular strain configuration consid- 
ered here would not be a necessary ingredient. Instead 
what is of utmost importance is the existence of a finite 
pseudo-magnetic field which could be achieved by any 
inhomogeneous strain configuration which generates the 
pseudo-Landau levels. 
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